Proceedings of the DIS'2004, Strbske Pleso, Slovakia 



FORWARD PHYSICS AT STAR 
STATUS OF ANALYSIS ON FORWARD AND MID RAPIDITY 
CORRELATION MEASUREMENTS IN P+P AND D+AU 

AKIO OGAWA FOR THE STAR COLLABORATION 

Brookhaven National Laboratory 
Upton NY 11973-5000, USA 
E-mail: akio@bnl.gov 

Measurements of the production of high energy tt" mesons at large pseudora- 
pidity, coincident with charged hadrons at mid-rapidity, for proton+proton and 
deuteron+gold collisions at -^/sjvjv = 200 GeV are reported. The p+p cross section 
for inclusive tt" production follows expectations from next-to-leading order pertur- 
bative QCD. Both the inclusive cross section and the di-hadron azimuthal correla- 
tions are consistent with a model where parton showers supplement leading-order 
pQCD. A suppression of the back-to-back azimuthal correlations was observed in 
d+Au, qualitatively consistent with the gluon saturation picture. 

It has been estabhshed from deep inelastic scattering (DIS) experiments that 
parton distribution functions (PDFs) are modified for nucleons bound in heavy 
nuclei. When the partons carry only a small fraction of the bound nucleon's mo- 
mentum (Bjorkcn x < 0.1), nuclear PDFs are found to be smaller than nucleon 
PDFsJL,. This phenomenon is known as shadowing. Establishing a quantitative 
understanding of nuclear PDFs is an essential step towards understanding the dy- 
namics of relativistic heavy ion collisions to address whether such collisions form a 
quark gluon plasma. Knowledge of nuclear gluon distribution functions is particu- 
larly important and is, at present, limited because of the limited range in resolution 
scale {Q^) spanned by the world data on nuclear DIS. Studies of proton-|-nucleus 
collisions at large center of mass energies {y^SNN=200 GeV) can provide constraints 
to the gluon density in heavy nuclei. 

In the pcrturbative QCD explanation of large-rapidity particle production, a 
large-x parton, typically a quark, scatters from a low- a; parton and then fragments 
into the observed particle(s). Forward charged particle production is found to be 
suppressed in d-f-Au collisions |2|, consistent with the expectation of gluon satura- 
tion 3 4 , indicating a different mechanism for particle production. Explanations of 
the suppression based on leading-twist pcrturbative QCD calculations employing a 
model of gluon shadowing have also been suggested |5j . Further tests of the possible 
role played by gluon saturation at RHIC energies could be provided by the study of 
particle correlations [H] • If the gluon density in the incident heavy ion is saturated, 
the large-x parton from the deuteron is expected to undergo multiple interactions 
through the dense gluon field resulting in multiple recoil partons instead of a single 
recoil parton. Hence, back-to-back di-hadron azimuthal correlations for d+Au col- 
lisions are expected to be smaller than those in p+p collisions. Coherent multiple 
parton scattering can also modify di-hadron azimuthal correlations 2l which may 
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Figure 1. Top view of the STAR detector. The scale is such that displacements transverse to the 
beam (x) are two times larger than displacements along the beam (z). 



be another approach to particle production in the presence of a saturated gluon dis- 
tribution. Leading twist calculations of particle production including conventional 
shadowing do not result in modifications of di-hadron azimuthal correlations jHj- 

An important question to address is whether fixed-order pQCD is appropri- 
ate to describe forward particle production in p-|-p collisions at •\/s=200 GeV. For 
< 62 GeV, next-to-leading order pQCD severely underpredicts measured 7r° 
cross sections At -^5=200 GeV and larger collision energies, there is quanti- 
tative agreement between NLO pQCD calculations and measured cross sections at 
mid-rapidity JOI- This agreement has been found to extend to 7r° production at 
(77) = 3.8 ^Jj. Further tests of the underlying dynamics responsible for forward 
particle production can be obtained from the study of particle correlations. In par- 
ticular, strong azimuthal correlations of hadron pairs are expected when particle 
production arises from 2^2 parton scattering, where in the hard scattering two 
initial-state partons scatter and result in two final-state partons. When a parton 
is scattered to fixed 77, the pseudorapidity distribution of the recoil parton is quite 
broad and is given by the convolution of PDFs with the angular distribution for 
partonic scattering. Extending pQCD beyond leading order introduces 2 — > 3 par- 
ton scattering that leads to a reduction of back-to-back azimuthal correlations of 
final-state hadron pairs. Similar effects are expected as the rapidity interval (Arj) 
between jet pairs increases |12| . 

This paper reports cross sections for forward inclusive tt'^ production for p+p 
collisions at y^=200 GeV. The azimuthal correlations between a forward tt" 
(('?7r)=4.0) and mid-rapidity charged hadrons were studied. These measurements 
for p-l-p collisions allow us to test our understanding of the particle production 
mechanism. In addition, exploratory studies with d-|-Au collisions at ^snn=200 
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GeV are reported, and azimuthal correlations of hadron pairs separated by large 
A77 were compared to those for p+p collisions. Forward tt*^ production in d+Au 
collisions refers to observation of the tt" in the direction of the incident deuteron. 

The Solenoidal Tracker At RHIC (STAR) is a multipurpose detector at 
Brookhaven National Laboratory. One of its principal components is a time projec- 
tion chamber (TPC) embedded in a 0.5T solenoidal magnetic field used for tracking 
of charged particles produced at <1.2. A forward tt" detector was installed in 
STAR as shown in Fig. ^ to allow for the detection of high energy tt" mesons at 
large rapidity. The FPD provides triggering and reconstruction of neutral pions 
produced with 3.3 < 7/ < 4.1. Data were collected over two years of RHIC opera- 
tions. In the 2002 run, p+p collisions were studied with a prototype FPD. Details 
about triggering, event reconstruction and normalizations are available in Ref. 
In the 2003 run, p-l-p collisions were studied with the FPD and exploratory mea- 
surements were also performed for d-|-Au collisions. 

The differential cross section for inclusive 7r° production for 30 < i?^ < 55 GeV 
at (77) = 3.8 was previously reported The event reconstruction and normaliza- 
tion methods were extended to allow measurement of the differential cross section 
at (77) = 3.3. The results are shown in Fig. [21in comparison to NLO pQCD calcula- 
tions evaluated at 77=3.3 and 3.8 using CTEQ6M parton distribution functions 
and equal renormalization and factorization scales of pt- The NLO pQCD calcula- 
tions are consistent with the data, in contrast to 7r° data at lower -^s P). The solid 
line was calculated using the "Kniehl-Kramer-Potter" (KKP) set of fragmentation 
functions ^1], while the dashed line uses the "Kretzer" set jTHj. The difference 
between the two reflects the uncertainty in the fragmentation functions at these 
kinematics. At the chosen scale, the KKP fragmentation functions tend to agree 
with the data better than Kretzer, consistent with what has been observed for 
midrapidity 7r° data at = 200 GeV [TH) . 

Correlations between a tt" produced at large rapidity with large Feynman x 
{xp > 0.25) and charged particles produced at I77I < 0.75 were also studied. The 
efhciency uncorrected average multiplicity of charged particles with |r7| <0.75 and 
Pt >0.2 GeV/c for events coincident with a forward tt" is 5.1 for p-l-p, and 14.4 
for d+Au. This multiplicity is approximately 10% larger than what is observed 
for minimum-bias events. The leading charged particle (LCP) analysis selects the 
mid-rapidity track with the highest px >0.5 GeV/c in each event. The azimuthal 
angle difference A0 — (p^o — 4>lcp is computed for each event. Distributions of A(/), 
normalized by the number of tt" observed at (77^)=4.0, for two tt" energy bins for 
p-l-p collisions are shown in Fig. |2| The left column of the figure shows simulations 
using PYTHIA 6.222 that account for detector resolution and reconstruction 
efficiency for both the forward tt*^ and the mid-rapidity charged particles. The right 
hand column are preliminary STAR data. 

The normalized At/) distributions were fitted by the sum of a constant and a 
Gaussian distribution centered at A0 = tt. Correlations near Ai^ = are not 
expected because of the large rapidity interval between the tt" and the LCP. The 
fitted parameters are the area under the Gaussian distribution (S), representing 
the azimuthally correlated coincidence probability; the azimuthally uncorrelated 
coincidence probability (B); and the width of the Gaussian (cr), having contributions 
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Figure 2. (Left) Inclusive vr" production cross section versus leading tt" energy (-Ejr) at average 
pseudorapidities {{rj)) 3.3 and 3.8. The average transverse momentum ((pt)) is correlated with 
En since the measurements were at a fixed r]. The inner error bars are statistical, and are smaller 
than the symbols for most data points. The outer error bars combine these with -dependent 
systematic errors added in quadrature. The curves are NLO pQCD calculations evaluated at r; = 
3.3 and 3.8 using different fragmentation functions as described in the text. (Right) Coincidence 
probability as function of azimuthal angle difference between tt" at forward rapidity ((»;)=4.0) 
and the leading charged particle at mid rapidity (|»7| < 0.75, PT > 0.5 GeV), in two different tt" 
energy bins for p+p collisions at -y^=200 GeV. The left column is for PYTHIA simulations and 
the right column is data. Fits to the data are described in the text. Error bars on the points 
represent only the statistical uncertainty. 



from transverse momentum in the hadronization of the jets and the momentum 
imbalance between the pair of jets {kr). The errors on the fitted parameters are 
based on the full error matrix. 

The PYTHIA simulation reproduces most features of the p+p data. The az- 
imuthally correlated di-hadron coincidence probability (S) can be identified in the 
PYTHIA simulation as arising from 2^2 parton scattering, resulting in a for- 
ward jet that fragments into the observed forward tt*^ and a midrapidity recoil jet 
that fragments into the observed charged hadrons. The width of the correlation 
may be underestimated by PYTHIA indicating that the average between the 
jet pair is too small. The azimuthally uncorrelated coincidence probability (B) 
primarily arises from 2^3 partonic processes, fully accounted for in NLO pQCD 
calculations, and approximately treated by initial- and final-state parton showers 
by PYTHIA. This interpretation is evident in the data when the azimuthal angle 
used in the correlation analysis is derived from the vector sum of the momenta of 
charged hadrons observed at midrapidity with pT >0.2 GeV/c with a further con- 
dition that the magnitude of the vector sum exceeds 0.5 GeV/c. The azimuthally 
correlated coincidence probability is found to be ^ 30% larger in such an analysis 
with a corresponding reduction in B. 

An exploratory data set for forward production was obtained for d+Au 



Forward Physics at STAR 



5 



71 + h~ correlations, simulation 
Vs = 200 GeV, Kr|^>| = 4.0, In J < 0.75 
p * p d + Au 



71 + h correlations, preliminary STAR data 

Vs = 200 GeV, Kr|^>| = 4.0, In J < 0.75 
p + p d + Au 




A^i = - 



S = 1 1 ± 1 % 




S = 2±1.1 % 




B = 50.4 ± 1 % 




: B = 88.1 ±1.6% 




- ag = 1.13 ±0.07 


0.2 


- Og = 0.54 ± 0.27 


<p,^>=1.1 GeV/c 








<Pt.lcp>=1" GeV/c 




0.1 






25<E<35 GeV 








2 4 6 




3 2 4 6 




Ai\> = 




*LCP 





Figure 3. Coincidence probability as function of azimuthal angle difference between tt" at forward 
rapidity ({tj) = 4.0) and leading charged particle at mid rapidity (|»7| < 0.75, PT > O.SGeV). (Left) 
HIJING 1.381 simulations for p+p and d+Au collisions, including detector effects. HIJING has 
substantially smaller fcy compared to PYTHIA 6.222, and it gives smaller a in p+p compared to 
PYTHIA (Fig. 0. (Right) STAR preliminary data for p+p and d+Au collisions. 



collisions. For this data, the calibration of the FPD is, at present, known to only 
10%. This prevents an accurate determination of the inclusive tt° cross section 
for d+Au collisions. Work is underway to improve upon this calibration. The 
slow variation of the p+p^ 7r° + correlations with 7r° energy suggests that 
similar observables in d+Au collisions should be relatively insensitive to calibration 
uncertainties. The normalized A(j) distribution for d+Au collisions, analyzed in the 
same way as the p+p data, is shown in comparison to p+p data in Fig. 13 

There is a large increase of BdAu relative to Bpp. This increase is also expected in 
default HIJING simulations 8 , as shown in Fig. |3| HIJING models d+Au collisions 
using PYTHIA for inelastic nucleon-nucleon interactions and the Glauber model 
to account for multiple collisions. The growth in B arises when a nucleon from the 
deuteron beam interacts with multiple nucleons from the Au beam. The HIJING 
simulation does not predict a significant difference between the width CdAu relative 
to CTpp. For the data, the fitted (XdAu is much smaller than Upp, most likely reflecting 
the inadequacy of the functional form used to represent A0dAu- The azimuthally 
correlated 7r° + coincidence probability is smaller in d+Au collisions than for 
p+p collisions. Part of the reduction is due to the fact that S is proportional to 
a and the fit results in a width for Ac/ijau that is likely to be unphysically small. 
HIJING simulations, that include a model of shadowing for nuclear PDFs, predict 
a significant azimuthally correlated tt'^ + coincidence probability. This is not 
observed in the preliminary STAR data. 

The largest systematic uncertainty in the Acf) distribution is expected to be the 
present understanding of the calibration of the calorimeter. Simulation studies of 
the calorimeter response suggest that an improved understanding of the calibra- 
tion can be obtained. The difference between Spp and SdAu is observed in several 
different methods of extracting an azimuthal angle from the charged particles ob- 
served with \ri\ <0.75. The impact of the functional form used to represent the 
Acf) distribution is under investigation. Complete assessment of systematic errors 
is underway. 

In summary, cross sections for the inclusive production of tt" mesons in p+p 
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collisions at -^5=200 GeV, at (77^)=3.3 and 3.8 are found to be consistent with 
NLO pQCD calculations, in contrast to the situation at smaller ^/s. The aziniuthal 
correlation between pairs of hadrons separated by large Arj is described by a leading 
order pQCD calculation including parton showers JHI for P+P collisions. Agree- 
ment of these calculations with the inclusive cross sections and di-hadron correla- 
tions strongly suggest that forward tt^ production arises from partonic scattering 
at this collision energy. This agreement means that forward particle production 
can be exploited as a probe of \ow-x parton densities in p-|-Au or d-|-Au collisions 
at y/SNN=200 GeV. Within the conventional theoretical framework, the lowest x 
values can be reached when both hadrons are detected in the forward direction 
|17|. The possible existence of gluon saturation in the Au nucleus may modify the 
correlation between the recoil jet rapidity and the initial-state Bjorken x value. 
Exploratory studies of forward 7r° production in d-|-Au collisions suggest that the 
azimuthally correlated component of hadron pairs separated by large rapidity in- 
tervals is suppressed relative to what is observed for p-|-p collisions and from model 
calculations. This is qualitatively consistent with the expectations that the particle 
production in a dense gluon medium differs from conventional leading-twist NLO 
pQCD expectations. More data for forward particle production and di-hadron cor- 
relations in d-|-Au collisions are required to reach a definitive conclusion about the 
possible existence of a saturated gluon state in the Au nucleus. A quantitative the- 
oretical understanding of the rapidity and pt distribution of di-hadron correlations 
would facilitate experimental tests of a possible color glass condensate. 
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